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Abstract. We discuss the role of the radiative B decays B — > X s ^y in our search for new physics 
. focusing on the issue of direct CP violation. We discuss in some detail a SM prediction for the CP 

asymmetries in inclusive b — > s/d transitions, namely | AS (B — > X s y) + As (B — > X^y) |~ 1 • 10~ 9 . 
Any measured value in significant deviation of this estimate would indicate new sources of CP 
violation beyond the CKM phase. 



in 



. Invited talks given by T.H. at the IPPP Workshop 'Phenomenology of Physics Beyond the SM', 

Durham (England), 6-11 May 2001, and at the International Workshop on QCD 'QCD@Work', 
ON , Martina Franca (Italy), 16-20 June 2001; to appear in the Proceedings. 

o 
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^ : INRODUCTION 

It is well-known that radiative B decays like B — > X s / d y play an important role in our 
search of new physics. These processes test the SM directly on the quantum level and, 
thus, are particularly sensitive to new physics (for a recent review see [1]). While the 
direct production of new (supersymmetric) particles is reserved for the planned hadronic 
machines such as the LHC at CERN, the indirect search of the B factories already implies 
significant restrictions for the parameter space of supersymmetric models and will thus 
lead to important clues for the direct search of supersymmetric particles. It is even 
possible that these rare processes lead to the first evidence of new physics by a significant 
deviation from the SM prediction, for example in the observables concerning direct CP 
violation, although it will then be difficult to identify in this way the new structures in 
detail. But also in the long run, after new physics has already been discovered, these 
decays will play an important role in analyzing in greater detail the underlying new 
dynamics. 

One of the main difficulties in examining the observables in B physics is the influ- 
ence of the strong interaction. As is well known, for matrix elements dominated by 
long-distance strong interactions there is no adequate quantitative solution available in 
quantum field theory. The resulting hadronic uncertainties restrict the opportunities in B 
physics significantly. If new physics does not show up in B physics through large devia- 
tions as recent experimental data indicates the focus on theoretically clean variables like 
inclusive radiative B decays is mandatory. 
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Within inclusive decays like B — > X s / d ythe long-distance strong interactions are less 
important and well under control due to the heavy mass expansion. The decay B — > X s y 
was first observed by the CLEO collaboration [2]; these measurements have been refined 
[3, 4] and confirmed by other experiments [5, 6] (see also [7]). The present world average 
using the present data from BELLE, CLEO and ALEPH is 

Q(B->X s y) = (3.22±0.40) x 1(T 4 . (1) 

The theoretical prediction of the Standard Model (SM) up to next-to-leading loga- 
rithmic (NLL) precision for the total decay rate of the B — > X s y mode [8] is well in 
agreement with the experimental data. The theoretical NLL prediction for the B — > X s y 
branching ratio [12] leads to 

<B(B^X s y) = (3.32 ±0.14 ±0.26) x 10~ 4 , (2) 

where the first error represents the uncertainty regarding the scale dependences, while 
the second error is the uncertainty due to the input parameters. In the second error the 
uncertainty due to the parameter m c /mb is dominant. 

This inclusive mode already allows for theoretically clean and rather strong con- 
straints on the parameter space of various extensions of the SM [9, 10, 1 1]. While many 
phenomenolgical analyses of the inclusive B — >■ X s / d y decays in supersymmetry are re- 
tricted to the assumption of minimal flavour violation including only CKM induced 
flavour change, a model-independent analysis also has to consider the generic new su- 
persymmetric sources of flavour change due to the mixing in the squark mass matrix. 
In [1 1] new model-independent bounds on supersymmetric flavour- violating parameters 
were derived from B — > X s y. The importance of interference effects for the bounds on 
the parameters in the squark mass matrices within the unconstrained MSSM is explic- 
itly demonstrated. In former analyses no correlations between the different sources of 
flavour violation were taken into account [10]. The new bounds are in general one order 
of magnitude weaker than the original bounds on the single off-diagonal element, which 
was derived in previous work [10] by neglecting any kind of interference effects. 

Recently, quark mass effects within the decay B — > X s y were further analysed [13], 
in particular the definitions of the quark masses m c and m b . The charm quark enters in 
specific NLL matrix elements (see Fig. 1) where the charm quark mass is dominantly 
off-shell. Therefore, the authors of [13] argue that the running charm mass should be 
chosen instead of the pole mass. The latter choice was used in all previous analyses 
[14, 8, 15, 16, 12]. 

m Poi 7m P° le m^(/i)/mP° le , //e[m c ,m,]. (3) 

Since these matrix elements start at NLL order only and, thus, the renormalization 
scheme for m c and m b is an NNLL issue, one should regard this choice as an educated 
guess of the NNLL corrections. The new choice is guided by the experience gained 
from many higher-order calculations in perturbation theory. Numerically, the shift from 

m? ole /m£ ole = 0.29 ± 0.02 to mf* (/j) / = 0.22 ± 0.04 is rather important and leads 
to a +11% shift of the central value of the B — > X s y branching ratio. With their new 
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FIGURE 1. a) Typical diagrams contributing to the matrix element of the operator O2 at the NLL level, 
b) Typical diagram contributing to the NLL anomalous dimension matrix. 

choice of the charm mass their theoretical prediction for the branching ratio is [13] 

<B(B^X s y) = (3.73 ±0.30) x 10~ 4 , (4) 

which induces a small difference between the theoretical and the experimental value. 

However, because the choice of the renormalization scheme for m c and is a NNLL 
effect, their important observation should be reinterpreted into a larger error bar in 
m^ s (jj) /m^ ole which includes also the value of irF c ole . A more conservative choice would 

then be m^ s (fj) /m^ ole = 0.22 ±0.07 which deletes the significance of the perceived 
discrepancy. 



DIRECT CP ASYMMETRY 

The CP violation in the B system will yield an important independent test of the SM 
description of CP violation (see [17]). In particular, detailed measurements of CP asym- 
metries in rare B decays will be possible in the near future. Theoretical predictions for 
the normalized CP asymmetries of the inclusive channels (see [18, 19, 20]) within the 
SM lead to 

^ B - ^ = r(^x, /rfY)+ r(^x f/jT ) (5) 

a CP (B -> X s y) « 0.6%, a CP (B -> B d y) w -16% (6) 

when the best-fit values for the CKM parameters [21] are used. An analysis for the 
leptonic counterparts is presented in [22]. The normalized CP asymmetries may also 
be calculated for exclusive decays; however, these results are model-dependent. An 
example of such a calculation may be found in [23]. 

CLEO has already presented a measurement of the CP asymmetry in inclusive b — » yy 
decays, yielding [24] 

a CP (B^X s y) = (-0.079 ±0.108 ±0.022) x (1.0 ±0.030), (7) 
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which indicates that very large effects are already excluded. 

Supersymmetric predictions for the CP asymmetries in B — > X s / d y depend strongly on 
what is assumed for the supersymmetry-breaking sector and are, thus, a rather model- 
dependent issue. The minimal supergravity model cannot account for large CP asymme- 
tries beyond 2% because of the constraints coming from the electron and neutron electric 
dipole moments [25]. However, more general models allow for larger asymmetries, of 
the order of 10% or even larger [26, 19]. Recent studies of the B — > gyrate asymmetry 
in specific models led to asymmetries between —40% and +40% [28] or —45% and 
+21% [27]. In general, CP asymmetries may lead to clean evidence for new physics by 
a significant deviation from the SM prediction. 

In [29] it was explicitly derived, that a bound on the combined asymmetries within 
the decays b — > sy and b — > dy, as well as their leptonic counterparts is possible which 
is more suitable for the experimental settings. It provides a stringent test, if the CKM 
matrix is indeed the only source of CP violation. Using U-spin, which is the SU (2) 
subgroup of flavour SU (3) relating the s and the d quark and which is already a well- 
known tool in the context of nonleptonic decays [30, 31], one derives relations between 
the CP asymmetries of the exclusive channels B — > K* y and B — > p~y and of the 
inclusive channels B — > X s y and B — > X d y. One should make use of the U-spin symmetry 
only with respect to the strong interactions. Moreover, within exclusive final states, the 
vector mesons like the U-spin doublet (K* , p~) are favoured as final states because 
these have masses much larger than the (current-quark) masses of any of the light quarks. 
Thus one expects, for the ground-state vector mesons, the U-spin symmetry to be quite 
accurate in spite of the nondegeneracy of m d and m s . Defining the rate asymmetries (not 
the normalized CP asymmetries) by 

Ar(fl -> V~y) = T(B -»• V~y) - T(B+ -> V+y) (8) 

one arrives at the following relation [29] : 

AT(B -> K*~y) + AT(B -> p~y) = b exc A exc (9) 

where the right-hand side is written as a product of a relative U-spin breaking b exc and 
a typical size A exc of the CP violating rate difference. In order to give an estimate of the 
right-hand side, one can use the model result from [23] for A exc , 

A exc = 2.5 x KT 7 T B . (10) 

The relative breaking b exc of U-spin can be estimated, e.g. from spectroscopy. This leads 
us to 

Mr* — m a 

=T777 ^ = 14%. (11) 

Certainly, other estimates are also possible, such as a comparison of / p and f K * . In this 
case one finds a very small U-spin breaking. Using the more conservative value for 
b exc , which is also compatible with sum rule calculations of form factors (see [32]), one 
arrives at the standard-model prediction for the difference of branching ratios 

| A$(B~ ^K*~y) +A®(5~ -> p~y) |~4x 10~ 8 (12) 
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Note that the right-hand side is model-dependent. 

Quite recently, the U-spin breaking effects were also estimated in the QCD factoriza- 
tion approach [33]. Within this approach, it was shown that the U-spin breaking effect 
essentially scales with the differences of the two form factors, (Fk* —F p ). Using the 
formfactors from the QCD sum rule calculation in [34] and maximizing the CP asym- 
metries by a specific choice of the CKM angle y, the authors of [33] obtain 

A$(B~ ^K*~y) +A<B(B- -> p~y) ~-3x 1(T 7 , (13) 

while for the separate asymmetries they obtain: 

A${B^K*i) = -7x 1(T 7 , A'B(B^py) =4x 1CT 7 . (14) 

This calculation explicitly shows the limitations of the relation (9) as a SM test. 

The issue is much more attractive in the inclusive modes. Due to the l/ra^ expansion 
for the inclusive process, the leading contribution is the free &-quark decay. In particular, 
there is no sensitivity to the spectator quark and thus one arrives at the following relation 
for the CP rate asymmetries [29] : 

Ar(5 - X s y) +AT(B - X d y) = b inc A inc . (15) 

In this framework one relies on parton-hadron duality; so one can actually compute 
the breaking of U-spin by keeping a nonvanishing strange quark mass. The typical size 
of bi nc can be roughly estimated to be of the order of | b[ nc | ~ m s /m b ~ 5 x 10 ; | A, nc | 
is again the average of the moduli of the two CP rate asymmetries. These have been 
calculated (for vanishing strange quark mass), e.g. in [18] and thus one arrives at 

| A<B (B -> X s y) + A<B(B^ X d y) |~ 1 x 1(T 9 . (16) 

Any measured value in significant deviation of (16) would be an indication of new 
sources of CP violation. Although only an estimate is given here, it should again be 
stressed that in the inclusive mode the right-hand side in (16) can be computed in a 
model-independent way with the help of the heavy mass expansion. 

Going beyond leading order in the I /nib expansion the first subleading corrections 
are of order 1 jm\ only. The 1 jm\ corrections are induced by the imaginary part of the 
forward scattering amplitude T(q): 

T(q) =i J d 4 x <B \ ToJ(x) O 7 (0) | B > exp(iqx) . (17) 

where only the magnetic operator O-j is taken into account. Using the operator product 
expansion for TO^(x) Oi(0) and heavy quark effective theory methods, the decay width 
T(B — ► X s y) reads [38, 39] (modulo higher terms in the l/m b expansion): 

(o l0l ) OiGlml , o o, , ( § NP \ 

5^ = ^1-^2 , (18) 
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(a) (b) 

FIGURE 2. a) Feynman diagram from which the operator O arises, b) Relevant cut-diagram for the 
(02, 0-])-interference. 



where X\ and X 2 are the parameters for kinetic energy and the chromomagnetic energy. 
Using X { = -0.5 GeV 2 and X 2 = 0.12GeV 2 , one gets 8 NP ~ -4%. Thus, the contribu- 
tions are small and cancel in the sum of the rate asymmetries - in the limit of U-spin 
symmetry. The U-spin breaking effects in this contribution also induce an overall factor 
m 2 /m^ in addition - as one can read off from the explicit results for the B — > X s l + l~ case 
[35]. 

There are also nonperturbative corrections which scale with 1/m 2 [36]. which are 
induced by the interference of the mangnetic O7 and the four-quark operator 0%. This 
effect is generated by the diagram in Fig. 2a (and by the one, not shown, where the gluon 
and the photon are interchanged); g is a soft gluon interacting with the charm quarks in 
the loop. Up to a characteristic Lorentz structure, this loop is given by the integral 

t 2 VTTT^ JrJr ■ (19) 

Jo jo m L c — k L g x\ 1 — x) — 2xyk g k^ 

As the gluon is soft, i.e. k 2 gl k g ky A^ C£) m^/2 m 2 , the integral can be expanded in k g . 
The (formally) leading operator, denoted by , is 

= ^f cb V*C 2 " j5)gsG vX be^d% a . (20) 

Then working out the cut diagram shown in Fig. 2b, one obtains the nonperturbative 

contribution T^_f^ y to the decay width, which is due to the (02,07) interference. 
Normalizing this contribution by the LL partonic width, one obtains 

^4pL--^^+0.03 . (21) 
r£, 7 9C 7 m 2 

This result corresponds to the leading term in an expansion in the parameter 
t = k g ky/2m 2 . The expansion parameter is approximately mfcAgc»/2m 2 ~ 0.3 (rather 
than Ag CZ) /m 2 ) and it is not a priori clear whether formally higher order terms in the m c 
expansion are numerically suppressed. However, the explicit expansion of the complete 
vertex function, corresponding to Fig. 2a, shows that higher order terms are indeed 
suppressed, because the corresponding coefficients are small (see i.e. [37]). 
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Moreover, the operator does not contain any information on the strange mass, thus, 
also in these contributions one finds the same overall suppression factor from the U-spin 
breaking. 

The corresponding long-distance contributions from up-quark loops are CKM sup- 
pressed in the B — > X s y case, but this does not hold in the B — > X^y case. Naively, 
one could expect that the corresponding contributions from up-quark loops scale with 
1/m^. However, following the approach of [37], one easily shows that the general ver- 
tex function cannot be expanded in the parameter t in that case. However, the expan- 
sion in inverse powers of t is reasonable. The leading term in this expansion scales like 
t~ l ~ m^/kgky and therefore cancels the factor l/m„ in the prefactor (see the analogous 
1/m^ factor in (20)) and one gets a suppression factor (Ag CZ) /m^) • (m^/k g ky) [37]. Thus, 
although the expansion in inverse powers in t induces nonlocal operators, one explicitly 
finds that the leading term scales with Agco /nib 1 - 

This argument improves the discussion in [29] where the argument was given that 
vector dominance calculations in [40] show that the long-distance contributions from the 
up-quark loops to the decay rates are found to be rather small. However, that argument 
in [29] does not allow any statement about the effects of the U-spin breaking in constrast 
to the one given here. 

Summing up, the analysis shows that the known nonperturbative contributions to (16) 
are under control and small and that this prediction provides a clean SM test, if generic 
new CP phases are present or not. Any significant deviation from the estimate (16) would 
be a strong hint to non-CKM contributions to CP violation. 

Finally, we emphasize that an analogous prediction for the leptonic inclusive B — > 
Xs/dl modes is also possible taking into account some specific cuts on the invariant 
dilepton spectrum. 
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